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In this lecture we discuss various properties of the phase factor of the fermion determinant for 
QCD at nonzero chemical potential. Its effect on physical observables is elucidated by compar- 
ing the phase diagram of QCD and phase quenched QCD and by illustrating the failure of the 
Banks-Casher formula with the example of one-dimensional QCD. The average phase factor and 
the distribution of the phase are calculated to one-loop order in chiral perturbation theory. In 
quantitative agreement with lattice QCD results, we find that the distribution is Gaussian with a 
width ~ \iT\fV (for m n <C T <C Aqcd)- Finally, we introduce, so-called teflon plated observables 
which can be calculated accurately by Monte Carlo even though the sign problem is severe. 
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1. Introduction 



Although the Euclidean formulation of QCD has been extremely successful in the evaluation 
of observables from first principles, at nonzero chemical potential, because of the phase of the 
fermion determinant, progress has been slow and limited to a small part of the QCD phase dia- 
gram (see for example [jl|] for a recent review). Motivated by the drastic effect of the phase of 
fermion determinant on physical observables, we study its properties to one-loop order in chiral 
perturbation theory. Since the average phase factor is a ratio of two partition functions it has to 
vanish exponentially with the volume when evaluated with respect to a positive definite measure. 
Therefore, a sensible strategy is to extract observables by extrapolation from small volumes. In 
the domain of validity of chiral perturbation theory it possible to derive the volume dependence of 
the average phase factor, and to determine the parameter range where reweighting methods [^, [|] 
work. The distribution of the phase, which is also calculated, might be relevant in applications of 
the density of states method [Q, ||, |6|]. It is less clear how methods that rely on analytical contin- 
uation in the chemical potential []7|, [8], ^, [H]] or methods that are based on the canonical partition 



function [11, 12] are affected by the sign problem, but each method has its own difficulties. 

With regards to the volume dependence of the QCD partition function at /I ^ 0, this can be 



investigated in detail in the microscopic domain of QCD [ 13 , |14|] where QCD is equivalent to chiral 



random matrix theory. In this domain, the cancellations that give rise to the discontinuity [ |15[ ] of 
the mass dependence of the chiral condensate also take place for finite size matrices Jl6|]. 

After the definition of the average phase factor in section 2, we discuss its effects on the 
phase diagram and on the chiral condensate in Id QCD in Section 3. Results obtained with chiral 
perturbation theory including the distribution of the phase are given in sections 4 and 5. 



2. Average Phase Factor 



i 



quark mass m 



At nonzero chemical potential, }X, the QCD Dirac operator is given by D + /i/o + m, where D 
is the anti-hermitian Dirac operator at zero chemical potential, and m is the quark mass. Because 
* /i/o is hermitian, this Dirac operator does not 

have any hermiticity properties and its eigen- 
values are scattered in the complex plane (See 
Fig. 1). Although the Dirac eigenvalues are 
still paired, the determinant is complex 
det(£> + /iy + m) = e w | det(D + jU% +m) \. 
The phase will be strongly fluctuating if the 
quark mass is inside the domain of eigenval- 
ues resulting in an average phase factor that 
vanishes in the thermodynamic limit. 
The average can be evaluated with respect to different partition functions. A particularly useful 
definition of the average is with respect to the phase quenched partition function, 




Fig. 1 . Quenched Dirac eigenvalues on a 4 3 
and a schematic Dirac spectrum (right). 
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which is the ratio of the QCD partition function and the phase quenched QCD partition function. 
Because free energies are extensive and the free energy of QCD and phase quenched QCD is 
generally different, the average phase factor vanishes exponentially with the volume. This explains 
the severity of the sign problem. 

In addition to the sign problem, there is an overlap problem. For any given field configuration, 
the chiral condensate is close to the result for the quenched or phase quenched partition function. 
The true value of the condensate is due to exceptional gauge field configurations. The deep reason 
for this is that determinants and the average phase factor depend exponentially on the volume. 



3. The effect of the phase factor 

The simplest way to analyze the effect of the phase factor is to compare the phase diagram 
of QCD and phase quenched QCD (denoted by IQCDI). First principle calculations for QCD have 
only been done at or near the temperature axis of the phase diagram in the jUT-plane. At zero 
temperature, a transition to a phase with nonzero baryon density takes place at /I m^/2>. Experi- 



mentally []18[], it has been shown that this transition is first order with the first order line ending in 
a critical endpoint (see Fig. 2). The remainder of the phase diagram is based on general arguments 
and model studies. At extremely high densities, one expects a superconducting phase whereas at in- 
termediate densities a confined phase with restored chiral symmetry is a likely possibility Jl9|, 




Fig. 2. The phase diagram of QCD in the temperature baryon chemical potential plane (left) and in the temperature isospin chemical 
potential plane (right). The data points in the right figure are from lattice QCD simulations [p"?|. 

On inspection of the lattice Dirac spectra at nonzero chemical potential (see Fig. 1) we notice 
that the eigenvalues are distributed more or less homogeneously inside a strip and that the strip 
has a sharp edge. This implies that we can expect a phase transition when the quark mass hits the 
edge of this strip located at A(ju). The critical chemical potential is therefore given by A(jU c ) = m 
[22, 23]. Because phase quenched QCD is QCD at nonzero isospin chemical potential, we know 
that at low temperatures a phase transition to a pion condensed phase occurs at /I = m n /2. Using 
the Gell-Mann-Oakes-Renner relation this results in a half-width of the strip of A(/x) = 2/x 2 F 2 /Z. 

The width Dirac spectrum for each gauge field configuration in full QCD is also given by 
2A(jii), but instead of a phase transition at pL = m n /2 there is a phase transition at /I = m^/3. 
In particular, this implies that the chiral condensate should be a smooth function of the chemical 
potential in the region with zero baryon density, and in the low temperature limit, it should not 
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depend at all on the chemical potential. How this can happen in spite of the fact that the Dirac 
spectrum is /I -dependent has been understood in detail for the microscopic domain of QCD [15 



16Q . Here we will illustrate this phenomenon for QCD in one-dimension [|24[]. 



3.1 Effect of the Phase in Id QCD 



The Banks-Casher relation [25] states that 



£ = lim lim 

m^OV^oc 



Tip [m) 
V 



(3-1) 



where p (A) is the density of Dirac eigenvalues on the imaginary axis and V is the volume of space- 
time. Although originally intended for an anti-Hermitian Dirac operator, this relation correctly 
gives a vanishing chiral condensate for phase quenched QCD at jj. / since there is no accumula- 
tion of eigenvalues on the imaginary axis. However, for full QCD at /I ^ 0, the chiral condensate 
has a discontinuity when the quark mass crosses the imaginary axis, but it does so without an accu- 
mulation of eigenvalues on the imaginary axis. The alternative mechanism that is at work was first 
discovered in random matrix theory [15, O]. Below we will illustrate this mechanism for QCD in 



one dimension which can also be viewed as a random matrix model. 











m 





Dirac spectrum of 1d QCD 

Fig. 3. Schematic plot of eigenvalues of the Dirac operator for lattice QCD in one dimension (left). The yellow dots denote the 
position of the eigenvalues for a single gauge field configuration, whereas the green ellipse shows the support of the spectrum in the 
thermodynamic limit. The chiral condensate for one flavor versus the quark mass is shown in the right figure. 

The partition function of lattice QCD in one dimension is given by 



UeU(N c ) 



dUdetD, 



(3.2) 



where the integral is over the Haar-measure of U (N c ). The Id lattice Dirac operator is the N xN c 
matrix with hopping matrix elements given by U exp(/i) and its inverse [ |26| ] (the number of lattice 
points is denoted by N). The chiral condensate for one flavor is given by 



E(jn) 



(3.3) 



Since the eigenvalues, Xt, are located on an ellipse in the complex plane (see Fig. 3) the determinant 



has a complex phase. For U(l) the chiral condensate can be evaluated analytically [24] with the 
result that is shown in the right figure of Fig. 3. The amazing phenomenon, also known as the 



"Silver Blaze Problem" [|27[], is that the chiral condensate is continuous when m crosses the ellipse 
of eigenvalues, but shows a discontinuity at m = 0, where there are no eigenvalues. This can happen 
because the chiral condensate is determined by exponentially large (in the number of lattice points) 
contributions which cancel to give a finite result for N — > °°. 
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4. Average Phase Factor in Chiral Perturbation Theory 



Because the average phase factor is the ratio of the QCD partition function and the phase 
quenched QCD partition function, it can be evaluated by means of chiral perturbation theory. The 
low-temperature limit of QCD is given by a gas of pions and does not depend on /I to one-loop 
order. Since the phase quenched partition function is QCD at non zero isospin chemical potential 
the pions are charged with respect to jU resulting in a jU -dependent one-loop result. For Nf = 2 there 
are two charged Goldstone modes n±. De- 
noting their one-loop contribution by 

1 



n 



Go(M)/2 



(4.1) 



TIT c = 0.76 

CPT 

Alllon 



p ^ml+p 2 + (p Q -2ilx) 2 
we obtain (the product is over all pions) 
<det 2 (D + m + J U7o)> ~ e~ VFm rj e G °C= )/ 2 , 





TIT = 0.90 
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The average phase factor is the ratio of 
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the two partition functions cf. (2.1) result- 
ing in the cancellation of the neutral Gold- 
stone bosons so that 



Fig. 4. Average phase factor in lattice QCD obtained 
from |0j (red curve) compared to one loop chiral per- 
turbation theory in a box equal to the size of the lattice 
(dashed curve) . 

In Fig. 4 where we compare one loop results for the average phase factor to lattice QCD results 
obtained by Taylor expansion in the chemical potential (there are no fitting parameters). The agree- 
ment is striking even for T/T c = 0.90. As expected chiral perturbation theory fails above T c . 



(e 2W ) 



pq 



„Go(ju=0)-Go(/t) 



(4.2) 



5. Probability Distribution of the Phase 

The density of the phase angle is defined by p(<j>) = (8(<p — Imlogdet(D + m + jU/o)))^. 
Notice that G (— oo 5 oo). If the average is over dynamical quarks, the phase density factorizes into 
an overall phase factor and the phase density of the phase quenched partition function (p pq (0)), 

(d^-d)e iN f e \det N f(D + m + ^Y )\) =e iN ^(d(^-d)\det N f(D + m + ^Y )\). (5.1) 

According to the Central Limit Theorem we expect that p pq (0) approaches a Gaussian distribution 



[28]. Notice however that observables are determined by correlations with the phase of the fermion 
determinant, and knowing the Gaussian distribution is only part of the story. 

The phase density can be evaluated using the replica trick and chiral perturbation theory 
Fourier transforming the 8 -function we obtain 



Pjv>(0) = (8{<j) -Imlogdet(D + m + ^7o)))i V/ = (£e 



in ( - Im log det (D+m+jl 70 ) > 



(5.2) 



The phase density therefore follows from the moments of the phase factor, 



(e 



2in9 \ 



1 / det" +yv /(D + m + Atyo) \ 



(5.3) 
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which is the ratio of two partition functions with a numerator that contains both bosonic and 
fermionic Goldstone particles. All 2n(n+Nf) charged Goldstone particles are fermions while 
the uncharged Goldstone particles are bosons. We thus find 



(e 



2ind\ 



,n(n+W / )[G (M=0)-Go(M)] 



(5.4) 



By Poisson resummation we obtain (AG = Go(ju) — Go(jU = 0)) 

Y e in$ e -n/2(n/2+N f )AG 



NjAG 



P(0) 



e*"f 



(5.5) 



A simple expression for the width of the Gaussian is 
obtained in a small fL expansion of the ChPT result. 
For the extensive quantity AG we find {m K <C T) 



AG 



V\x 2 T 2 . 



(5.6) 



Since the width of the Gaussian ~ yAG this result 
is in quantitative agreement with the lattice results 
by Ejiri |t§] (See Fig. 5). 




9=(JV/4)Im[ln(det A/)] 

Fig. 5. Phase density in lattice QCD g]. 
The average phase factor can also be evaluated in a weak coupling expansion valid at high temper- 
atures 



31]. The small-^t limit of the first nonvanishing contribution is given by [ |29| ] 

F QCD (n,T)-F ]QCD] (n,T) 



Vafn 2 T 2 , 



(5.7) 



with V the space-time volume. This implies that the sign problem remains severe in the thermody- 
namic limit of the weak coupling result. Notice that expectation values can be calculated accurately 
despite the sign problem if the correlations between the observable and the phase factor are small. 
Such teflon plated observables should obey the condition (here Nf = 2) 

< 1 



g2 ' e )|QCD| 



| QCD 



(e 2 



| QCD | 



This condition only makes sense if 



QCD ( e2 ' e ) QCD 

qcd| 7^ as can be seen by rewriting the l.h.s. as 

|(^)qcd -(^)|qcd|I 

I QCD I 



(5.8) 



(5.9) 



For example, when G is the baryon density this ratio is 0(a 2 ) whereas (^)qcd is 0(a®) 



6. Conclusions 



The phase factor of the fermion determinant completely alters the properties and the physics 
of the QCD partition function at nonzero chemical potential. We have illustrated this by the rela- 
tion between the chiral condensate and the Dirac spectrum and the difference between the phase 
diagram of QCD and phase quenched QCD. Using one-loop chiral perturbation theory, we have de- 
termined both the magnitude of the average phase factor and the statistical distribution of the phase 
which turns out be a Gaussian with width ~ [iTy/V. Finally, we introduced a class of observables 
(teflon plated observables) that can be evaluated accurately despite the sign problem. An example 
of such observable is the baryon density at high temperature and small chemical potential. 
Acknowledgments: This work was supported in part by U.S. DOE Grant No. DE-FAG-88ER40388. 
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